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We have updated our previous investigation of the production of muons and neutrinos in cosmic ray interactions 
with the atmosphere, taking account of recent results from the ep collider HERA in our QCD-based model for 
hadronic interactions. Qualitatively, our previous results remain unmodified: our predictions for the conventional 
muon and neutrino fluxes agree with earlier calculations, whereas the charm particle treatment we use gives 
significantly lower prompt fluxes compared to earlier estimates. This implies better prospects for detecting very 
high energy neutrinos from cosmic sources. 



The flux of muons and neutrinos at the Earth 
has an important contribution from decays of par- 
ticles produced through the interaction of cosmic 
rays in the atmosphere. This has an interest in its 
own right, since it reflects primary interactions at 
energies that can by far exceed the highest avail- 
able accelerator energies. It is also a background 
in studies of neutrinos from cosmic sources as at- 
tempted in large neutrino telescopes. 

We have updated our detailed study of muon 
and neutrino production in cosmic ray interac- 
tions with nuclei in the atmosphere by taking 
into account the information on nucleon structure 
functions from recent experiments at the ep col- 
lider HERA [^IJ] and by refining our simulation 
of particle cascades in the atmosphere with a bet- 
ter treatment of secondary particles. A complete 
description is presented in ref . § . 

Our investigations agree with earlier studies 
on the fluxes of atmospheric muons and neu- 
trinos coming from decays of tt and K mesons. 
On the other hand, we flnd quite low fluxes for 
prompt muons and neutrinos, which arise through 
semi-leptonic decays of hadrons containing heavy 
quarks (most notably charm) . Other estimates of 
these prompt fluxes [^-|ll| are higher than ours, 
and vary by few orders of magnitude. These dif- 
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ferences come from extrapolating charm produc- 
tion data from accelerator energies to the orders- 
of-magnitude higher energies of the relevant cos- 
mic ray collisions. Current data from surface and 
underground detectors attempting to measure the 
flux of prompt muons and neutrinos (see e.g. p2[) 
are still too discrepant to discriminate between 
the different models for charm production. 

The main contribution of our study is in us- 
ing proper charm production data and a sound 
physical model based on QCD. First, we use re- 
cent charm cross section measurements that form 
a consistent set of data, but disagree with some 
of the early measurements that were substantially 
higher. Secondly, we apply state-of-the-art mod- 
els to simulate charm particle production in high 
energy hadron-hadron interactions. 

We have obtained the atmospheric muon and 
neutrino fluxes with two different methods: via a 
Monte Carlo simulation of the hadronic cascade 
and via approximate analytical expressions with 
energy-dependent Z-moments. We were satisfied 
that the two methods, which are conceptually 
rather different, gave quite similar results. Dif- 
ferences were typically less than 20%, below the 
uncertainty in our charm calculation and quite 
acceptable in this context. 

Very different prompt fiuxes are predicted by 
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Figure 1. Our prompt and conventional muon and neutrino fluxes (solid lines) compared with those in 
other models for charm production. 



other models for charm production, which differ 
both in the magnitude and energy dependence of 
the total charm production cross section and in 
the distribution in longitudinal momentum frac- 
tion of the charmed particles. 

Volkova et al. applying the so-called 'quark- 
gluon string model' p^ , obtained the curves la- 
beled V in fig. |l|. CastagnoH et al. obtained the 
result marked C. Both models use a parametrized 
energy dependent charm cross section (curves C 
and V respectively in fig. ^) normalized to early 
experimental data which are substantially above 
later measurements. A more recent calculation 
based on the QGSM by Gonzalez-Garcia et al. |^ 
gives fluxes that are comparable to curve C in 

Curves marked Zl and Z2 in figs. |^ and ^ are 
from Zas et al. Curves Zl illustrate an ex- 
treme model where the charm cross section is 
simply taken as 10% of the total inelastic cross 
section (which is curve cr*°* in fig. H). This is sub- 
stantially higher than all charm data. Curves Z2 
correspond to charm quark production calculated 
with leading order perturbative QCD matrix ele- 
ments using relatively hard parton distributions. 

All the previous models except Z2 assume 
Feynman scaling for the charm energy spectra. 
Bugaev et al. ||l^ considered Feynman scaling 



violations through a phcnomcnological equation. 
Their calculation resulted in overall prompt fluxes 
slightly lower than Volkova's (curve V). Their 
fluxes are higher in the non-scaling case than in 
the scaling case (opposite to our results). 

Our model (solid lines in fig. |l| and curve MC in 
fig. 1^) uses leading order perturbative QCD ma- 
trix elements supplemented by a correction K = 2 
for next-to- leading order processes. It incorpo- 
rates a strong breaking of Feynman scaling, aris- 
ing from the dominance of perturbative charm 
quark production close to the charm threshold. 
For the parton densities in the nucleon we have 
adopted the MRS G parametrization , which 
uses essentially all relevant experimental data, 
from deep inelastic scattering experiments to re- 
cent results from the ep collider HERA 
Actually, since a naive extrapolation of the G 
parametrization below the measured region in x 
at rather small "^c leads to a presumably 

unphysically large charm production cross sec- 
tion, we have implemented a flatter dependence 
as x^O like x'' with e ~ 0.08 (cfr. |l|). Curves 
labeled G and Dq in fig. || show the charm produc- 
tion cross section one would obtain from the (un- 
fiattened) G parametrization and from the now- 
unacceptable MRS Do parametrization 

Some of the previous models give charm pro- 
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Figure 2. Energy dependence of charm produc- 
tion cross section in pp (pp) collisions. Refer- 
ences to experimental points are given in ref . . 
Curves are described in the text. 

duction cross sections higher than recent data 
(see fig. H) and some assume Feynman scaling 
for the longitudinal momentum distributions. In 
some cases one may have been mislead in the 
construction of the models by the early charm 
measurements that turned out to be substantially 
higher than the measurements done later. Our 
model instead gives a fair description of mea- 
sured charm production cross sections and applies 
well-motivated charm particle momentum distri- 
butions with significant Feynman scaling viola- 
tions. 

With respect to other models for charm pro- 
duction, we predict substantially lower muon and 
neutrino atmospheric fluxes from decays of charm 
mesons. In particular, we predict an interest- 
ingly low prompt neutrino background to search- 
ing high energy neutrino cosmic sources in large 
scale neutrino telescopes. 
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